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We present preliminary results of a study of A5 = 2 matrix elements originating from physics 
beyond the Standard Model. Using 2 + 1 flavour Domain Wall Fermions we obtain the non- 
perturbative renormalisation (mixing) matrix in the RI scheme. We also discuss plans for the 
chiral extrapolation of the renormalised matrix elements in a partially quenched set-up. 
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1. Motivation 



In the Standard Model (SM) there is only one four-quark operator contributing to the analysis 
of Kq — Kq mixing. It originates from box diagrams with W bosons, Figure [l[ In theories beyond 
the SM there are additional box diagrams. Figure ^ shows an example from the mass insertion 
approximation of a supersymmetric extension of the SM. 





Figure 1: Box diagram for Kq-Kq mixing in the 
Standai-d Model. 



Figure 2: Example for a SUSY box diagram 
with gluinos and squarks in the mass insertion 
approximation. 



The matrix element of the SM operator 

Qi=rynPLd"ffPLd\ (1.1) 

is parametrised by the bag parameter Bk- Q\ renormalises multiplicatively for lattice discreti- 
sations with chiral symmetry. The most advanced computations of Bk have exploited the chiral 
symmetry of Domain Wall Fermions (DWF) [1] or overlap fermions [2]. The good agreement be- 
tween the measured value for |e| and the SM with lattice Bk input puts a constraint on theories 
beyond the Standard Model (BSM) which allow for additional classes of box diagrams. The lattice 
can contribute here as well by providing hadronic matrix elements for a basis of four-quark opera- 
tors which then can be used in phenomenological studies. Again lattice fermions with good chiral 
symmetry are essential since renormalisation becomes continuum-like. 
Besides the SM operator Qi (|0|) we consider the operators, 



Q2 = rY^PLd"ty^'PRd\ (1.2) 

Q,=rPLd"s^PRd', (1.3) 

Q4=rPLd"s'PLd', (1.4) 

Q5=ra^yPLd"s'a^'PLd''. (1.5) 



This basis is related to the supersymmetric basis used in several other studies by a Fierz transfor- 
mation. The anomalous dimensions for these operators are known to next-to-leading order [3,4]. 
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The mixing among the BSM operators is Umited to two 2x2 blocks, 
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2. Non-perturbative Renormalisation 

To subtract the logarithmic divergence in the matrix elements of the operators Qi we employ 
the RI scheme [5] at a scale fi. We consider the full 5x5 matrix, 

Qf = Zij{^)Qj, (2.1) 

despite the reduced mixing of the continuum theory. Since our lattice action is chirally symmetric 
up to a small violation of O(mres) we expect the elements of Z which vanish in the presence of 
chiral symmetry to be small. Our approach enables us to check this assumption and to compare the 
renormalised results with the full matrix and the block diagonal one. 

We compute the amputated four-point vertex functions F and project them onto the relevant 
Dirac structure. The resulting matrix A yields the renormalisation matrix, 

i^ZlM) = Atree -A-'l/ = m')- (2-2) 

The wave function renormalisation Zq is eliminated using the (local) axial current and its renor- 
malisation constant Za which can be determined independently from the axial Ward identity for 
Domain Wall Fermions [6], 

^Z = Atree-A-VAA. (2.3) 

The standard RI renormalisation condition is defined at an unfortunate kinematic point. The 
momenta of all four involved quarks are the same and there is no momentum flowing out of the 
vertex. This so called exceptional momentum configuration leads to additional spontaneous chiral 
symmetry breaking from subgraphs which are only suppressed as a small inverse power of momen- 
tum. This has been seen in the splitting between and Ay, the amputated vertex functions for the 
axial and vector current [6]. The effect on the renormalisation factors for the operators Q2 to Q5 
is a ^ divergence. There are at least two suggestions in the literature how to deal with this pion 
pole. One can try to directly fit the pole in A;y [7] or one can form ratios of the type [8] 

miAij{mi , p^) - m2Aij{m2,p^) ^^4) 

m\—m2 ' 

where rui and ma are two different quark masses. The ratio in the chiral limit corresponds to the 
subtracted A^- [9]. 

We have implemented a slightly different renormahsation condition which directly gives the 
amputated Green's function without the pion pole. In this non-exceptional case one still has pf = 
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/i^ for all legs of the vertex, but there are two pairs of momenta which differ such that the sum of 
momenta at the vertex is also /i^. Now chirality breaking subgraphs are suppressed by the relatively 
high scale. 

In Figure |3| we compare the two kinematic set-ups for the row of the matrix A which determines 
the mixing for Q2. The main difference is the suppression of mixing with the operators ^4 and Q$ 
which is caused by the mentioned chirality breaking subgraphs. This clearly shows the advantages 
of the non-exceptional momentum data. The main obstacle for the use of this new approach is the 
lack of perturbative results for the matching to perturbative schemes like MS. A first result has been 
presented at this conference for the case of Z,„ [10]. Such infrared problems not only appear in the 
non-perturbtative data, but also seem to influence the rate of convergence of continuum perturbation 
theory, which is very promising for the RI-MOM approach. 
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Figure 3: Elements of the matrix A/A^ associated with the mixing for operator Q2 for a range of with 
exceptional momenta (left) and non-exceptional momenta (right). The statistical errors are smaller than the 
symbols. The chiral symmetry of DWF relating Q\ and Q2, Q3 and Qa, is manifest in the non-exceptional 
case, but obscured by spontaneous chiral symmetry breaking at low in the conventional approach. 

Recently we have adopted the use of gauge-fixed momentum sources which project on a single 
momentum at the source. The vertex then has to be put at the sink and a full volume average 
is possible. These sources allow us to reach a better statistical accuracy at much reduced cost 
especially on large volumes. Even though each momentum requires new propagator inversions the 
new method scales much better with the lattice volume and has to be the method of choice for our 
new large lattices. 

3. Results for Bare Matrix Elements 



The preliminary results presented in this section are part of the ongoing project with 32^ x 64 x 
32 ensembles of A/ = 2 + 1 Domain Wall Fermions with the Iwasaki gauge action at j3 = 2.25 [11]. 

The SM matrix element is commonly normalised by its value in the vacuum saturation approx- 
imation. Using the same normalisation for the other four operators leads to a divergent chiral limit 
for the resulting bag parameters. This divergence can be moved into the normalisation, but this 
requires the knowledge of renormalised quark masses for the phenomenological use of the results. 
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Therefore alternative normalisations like are preferable. But for the purpose of this write-up we 
stick with the Bk normalisation since all results are preliminary, 

{K\Qi\K)=NmlFlBi, /=1,...,5, (3.1) 
8 4 5 

Ni = -,--RaR,-R,-4R, (3.2) 

^=(^X- (3-3) 

Gauge-fixed wall sources have been shown to be a very efficient for B/f [1, 12]. We reach the 
same level of statistical uncertainty when applying this approach to the whole operator basis. In 
Figure ^ we give an impression of the quality of the data for the four non-SM operators. Each plot 
shows the fully dynamical data (myai = /Wsea) for both light and strange quark mass. With around 
100 configurations the statistical errors are of the order of 1% at the lowest dynamical mass for the 
bare B parameters. 



Our plan for the chiral extrapolation of this data set is again along the lines of the existing 
work on Bk- We use an approach based of SU (2) Chiral Perturbation Theory (ChPT). In this setup 
the leading order Low Energy Constants F and Bq are fixed from fits in the pion sector. The kaon 
is treated as a heavy meson [13]. The needed partially quenched formulae can be derived from 
results in Heavy Meson ChPT [14]. In Figure ^ we show the light quark mass dependence for 
BSM operators. 

4. Summary 

We have given a status report on a computation of AS' = 2 matrix elements for a complete op- 
erator basis. The aim of this study is to compute the renormalised matrix elements with dynamical 
fermions to have an impact on phenomenological studies. 

The use of Domain Wall Fermions with very small chiral symmetry breaking is essential for 
the continuum-like renormalisation of the operators. We use the RI scheme with gauge-fixed mo- 
mentum sources which greatly improves the statistical accuracy of the renormalisation constants. 
We are exploring a variant of the RI scheme with a non-exceptional momentum configuration to 
reduce effects from chiral symmetry breaking at low momentum. 
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Figure 4: Plateau plots for the ratios B2-B5 at the unitary quark masses of all three ensembles without 



normalisation factor A^,- (3.2) 



We have computed the bare matrix elements on the 32^ x 64 x 16 ensembles with small statis- 
tical errors for six valence masses. We plan to extrapolate this data using SU(2) Chiral Perturbation 
Theory in the same way as for Bk- 

We also intend to extend this work to the coarser 24^ lattices which will allow us to quantify 
the size of 0{a^) lattice artefacts. 
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Figure 5: Light quark mass dependence for the bag parameters B2-B5. The different colours indicate the 
three sea quark masses. In each plot there is data for two strange quark masses, 0.025 and 0.03. 
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